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Abstract 



, The well known domain wall type solutions are nowadays of great phys- 

■ ical interest in classical field theory. These solutions can mostly be found 

' only approximately. Recently the Hilbert-Chapman-Enskog method was 

< ^* ) succesfully applied to obtain this type solutions in 4> 4 theory. The goal 

of the present paper is to verify these perturbative results by numerical 
computations. 



Introduction 



Oh. During the last few years domain walls, strings and vortices have been of 



great physical interest in classical and quantum field theory. These are the 
solutions of highly non-linear partial differential equations with the energy con- 
centrated in a small region of the whole space-time. They play an important 
, role in cosmology, condensed matter and elementary particles theory. Usually 

t-j ' the exact form of these solutions can not be obtained as the appropriate equa- 

tions are too difficult to integrate. There are many approximate methods which 
can help to solve these equations but they usually do not provide the quanti- 
tative estimation of errors. Thus they should be compared with the numerical 
solutions obtained with help of the well tested algoritms. Numerical methods 
can also verify many assumptions which must be made in perturbative calcu- 
lations. In the present paper we consider the domain wall type solutions in </> 4 
theory. We compare our numerical results with the approximate solutions ob- 
tained in [Q . These approximate solutions were found with help of the so-called 
Hilbert-Chapman-Enskog method [0. It is the singular perturbative scheme || 
as the small perturbation parameter appears in the term with the highest order 
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derivatives. Thus from the mathematical point of view the problem is very sub- 
tle and its numerical verification is especially valuable. Let us stress that the 
paper is the part of the larger program ||, [|, [5] devoted studying different 
approximate methods in classical field theory. Many technical problems have 
been examined in this program so far and that is the second reason why we 
extend it on numerical analysis. 

Approximate solutions 

In this section we shortly describe the approximated domain wall type solu- 
tions in theory proposed in ^] . The lagrangian density and the corresponding 
field equation have the form 

/ m 2 \ 

d^d^- 2A ( $ 2 - — J $ = 0. (2) 

Here (j]^) = diag(— 1, 1, 1, 1) is the metric tensor and A , M are some positive 
constants. There exist two classical vacuum solutions in this model, namely 
$ = ±M/2\/A. A domain wall type solution smoothly interpolates between 
these two vacua. It is convenient to introduce in Minkowski space the special 
coordinate system (u°, u 1 , u 2 , s) co- moving with the domain wall |(|. The 
three coordinates u°, u 1 , u 2 describe the world volume of the membrane attached 
to the domain wall and £ — 2s /M is the distance from this three-dimensional 
hypersurface ( for details see ). Thus 

x» = Y> l (u)+€n fl (u), (3) 

where the functions Y^(u) describe the above mentioned hypersuface and 

is the normalised four- vector orthogonal to it. In these new coordinates the 

equation of motion takes the form 



M 2 ^Gdu a \ du b J 2ds 2 2h ds ds 

where 



$M = ^K,.) 1 (5) 
G ab = N ac g cd N db , N ac = g ac ~ £K ac , (6) 
r -G = V^gh.g = det[g ab ],G = det[G afc ]. (7) 
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The induced metric g a b and the external curvature K ac are defined by 



Kab =n lx d a d b Y»,i 



9ab = daY^dbY^. 



(8) 



The Hilbert-Chapman-Enskog method has been applied in |jj to obtain the per- 
turbative solution of the equation (4) . According to this method the following 
form is assumed for the approximate solution of the equation (4) 



In the present paper this approximation is verified numerically to give very 
precise results. However in general the problem is a little bit more complicated 
as the equation (4) involves the small perturbation parameter in the term with 
the highest order derivative. Usually in such case one can expect the solution to 
involve some extra non-analitic terms which of course are absent in the formula 
(9). This problem is discussed mainly from the physical point of view in JjJ. 
The non-analitic terms can really exist in the general solution but they are only 
radiative corrections to the base domain-wall configuration. These terms can 
be completely removed from the solution by the appropriate choice of the initial 
conditions. This is the case we are especially interested in as we would like 
to extract the pure domain-wall type solution without radiation. In practice 
it is very difficult to find the desired initial conditions. However even if the 
radiative corrections are present they quite rapidly decrease with time as the 
whole solution tends to the stable domain-wall type configuration. Thus these 
possible extra terms can practically be neglected. This is the reason why the 
Hilbert-Enskog-Chapman method could be safely applied in jjj. We will return 
to this problem in our conclusions. Assuming that for u° = the surface = 
coincides with the membrane s — the result looks as follows 




(9) 



4> = tanh s + 



1 C 



M cosh 2 s 





+ / 3 (s)4 2) (*) +-, 



(10) 



where 



(2) 



- sinh(2s) + - tanh s + - 



s 



(11) 



'2 



8 cosh 2 s ' 
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(12) 




(13) 
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f 3 (s) = - / dx—j- (14) 
Jo cosh x 

and the function C(u a ) satisfies the equation 

a^c + (j- - l) K b a KSC + ^C 3 = (y - i) ^X c , (15) 

with initial conditions 

BC 

C(u° = 0) = —(u° = 0) = 0. (16) 

Here the symbol cK 3 ) denotes the three-dimensional d'Alembertian on the 
membrane surface. The function C(u a ) determines the position of the surface 
where the field <j) vanishes. This position is given by the formula 

s = -^ + 0(M- 3 ). (17) 



Numerical results 

We have examined in detail two special cases of the general result (10), 
namely the cylindrical and the spherical domain walls. In the cylindrical case 
the function C(u a ) and the invariant K^K^K^ vanish. The other invariant 
K^K^ takes the simple form 



^X a = 2-f, (18) 



where 



r(tt°) = ro cos — (19) 

and ro is the initial position of the domain wall. Let us stress that in this 
case the position of the surface where the field <p vanishes coincides with the 
membrane attached to the domain wall, at least up to the second order in 1/M. 
The equation of motion (4) for the cylidrical domain wall reads as 




ro J \r r 

M 2 



1) = 0. (20) 
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The spherical case is a little bit more complicated because both invariants K^K h a 
and K^K^Ki are different from zero 



,,.4 
I n 



K b a K§ = 6-f , (21) 
KKlKl = S, (22) 



r 



r(u°)=r cn(V2 — ,l) (23). 



where 



r ' 2, 

As a consequence the function C(u°) does not vanish and should be obtained 
numerically. Thus in this case the distance between the above mentioned sur- 
faces is different from zero. The appropriate equation of motion for the spherical 
domain wall takes the form 




= 0. (24) 

We have solved the equations (20) and (24) numerically for —0.5 < £ < 0.5 and 
< u° < 1 with the initial data 0| u o =o , d(f>/du a \ u o = Q computed from formula 
(10) with KgK^ given by (18) or (21) respectively. Let us note that the maximal 
value of u° can not be too big because of possible singularities in coordinate 
system Q . We have introduced a grid of 8000 points and the partial derivatives 
with respect to s have been replaced with 

d(f) _ <t> n +l ~ 4>n-l d4> 2 _ 4>n+l - 20n + <f>n-l ^5) 



ds 2h d 2 s h 2 

The system of ordinary differential equations obtained this way has been inte- 
grated by the Runge-Kutta method ||. The numerical computations in both 
cases e.g. cylindrical and spherical have been performed for 7' = 2.5 and three 
different values of the parameter M = 10,20,40. We have compared the nu- 
merical functions (j> n um{u° — 1.0) and d<fi num / 'du° '(u° = 1.0) with their ap- 
proximate counterparts <p a p{u° — 1.0) and d<p ap / dvP(v? = 1.0) obtained from 
formula (10). The motion of the surface </> — has been also computed and 
compared with the approximate result (17). Figures (l)-(4) show the differ- 
ences A0 = 4>num — 4>ap and A(f> u o — d4> num /dvP — d(f> ap /du in the case of 
the cylindrical domain wall. We do not present the functions (f> num themselves 
because the interesting higher order corrections are very small compared to the 
u° independent term </>(°) = tanh(s). Instead the difference 4>num ~ tanh(s) is 
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presented in Fig(5) for M = 20. The corresponding derivative d<j)/du° is shown 
in Fig(6). The motion of the core 4> n um — is plotted in figures (13) - (15). 
Let us recall that in this case C = and the surface 4> ap = coincides with the 
membrane s = 0. The differences A<j> and A.<f> u a in the spherical case are pre- 
sented in figures (7) - (10). Figures (11) and (12) are the spherical counterparts 
of (5) and (6). The motion of the core <fi num = and the correspondig results 
computed from the formula (17) for <j) ap — are plotted in figures (16) - (18). 

Conclusions 

The aim of our paper was to verify the approximate solutions proposed in |Q 
by the numerical computations. Our investigations prove that the coincidence 
of the numerical results and the approximate ones is quite good. The existing 
differences could be affected by radiative corrections absent in the perturbative 
scheme. In numerical calculations these corrections can not be completely re- 
moved. However our results show that these possible terms are small and can 
be neglected in our considerations as the perturbative solution is only known up 
to the second order term in l/M. It would be interesting to compare these ra- 
diative corrections with the higher order terms of the perturbative scheme. This 
is the argument for further studies of these important questions. As it should 
be expected the differences 4> n um — 4>ap and d(/) nurn /du — d(f>ap/du Q decrease 
rather rapidly when the models with thinner walls (bigger M) are considered. 
Let us also stress that the numerical and the approximated motion of the do- 
main wall almost coincide for our biggest value of the parameter M. Thus the 
approximation proposed in |l]] seems to be an efficient tool to examine important 
details of some peculiar solutions in classical field theory like domain walls or 
vortices. In general our paper proves that the progress in the practical dealing 
with non-linear partial differential equations can be achieved when one combines 
different clever approaches such as the Hilbert-Enskog-Chapman perturbative 
method and the suitable choice of the coordinate system. 
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Fig. 1. The differences A<j)(u° = 1.0) for the cylindrical domain wall, M = 10 
and M = 20. 

Fig. 2. The differences A^> u o(u° = 1.0) for the cylindrical domain wall, 
M = 10 and M = 20. 

Fig. 3. The differences A0(u° = 1.0) for the cylindrical domain wall, M = 20 
and M = 40. 

Fig. 4. The differences A</> u o(u° = 1.0) for the cylindrical domain wall, 
M = 20 and M = 40. 

Fig. 5. The difference ^ mm («° = 1.0) — tanh(s) for the cylindrical domain 
wall, M = 40. 

Fig. 6. The derivative d4>/du°(u° = 1.0) for the cylindrical domain wall 
M = 40. 

Fig. 7. The differences A</>(u° = 1.0) for the spherical domain wall, M = 10 
and M = 20. 

Fig. 8. The differences A(f> u a(u° = 1.0) for the spherical domain wall, 
M = 10 and M = 20. 

Fig. 9. The differences A(f>(u a = 1.0) for the spherical domain wall, M = 20 
and M = 40. 

Fig. 10. The differences A</> u o(it° = 1.0) for the spherical domain wall, 
M = 20 and M = 40. 

Fig. 11. The difference 4> n um(u° = 1.0) — tanh(s) for the spherical domain 
wall, M = 40. 

Fig. 12. The derivative d(j>/du (u — 1.0) for the spherical domain wall 
M = 40. 

Fig. 13. The motion of the surface 4> n um = for the cylindrical domain 
wall, M=10. 

Fig. 14. The motion of the surface 4> nU m = for the cylindrical domain 
wall, M=20. 

Fig. 15. The motion of the surface (j)num = for the cylindrical domain wall 
M=40. 

Fig. 16. The motion of the surfaces 4> nU m — and (f> ap — for the spherical 
domain wall, M=10. 

Fig. 17. The motion of the surfaces </> TlUTO = and (f> ap — for the spherical 
domain wall, M=20. 

Fig. 18. The motion of the surfaces 4> nU m — and <j) ap — for the spherical 
domain wall, M=40. 



